Fitting Martingales To Given Marginals 

George Lowther 
Abstract 

We consider the problem of finding a real valued martingale fitting 
specified marginal distributions. For this to be possible, the marginals 
must be increasing in the convex order and have constant mean. We 
show that, under the extra condition that they are weakly continuous, 
the marginals can always be fitted in a unique way by a martingale 
which lies in a particular class of strong Markov processes. 

It is also shown that the map that this gives from the sets of 
marginal distributions to the martingale measures is continuous. Fur- 
thermore, we prove that it is the unique continuous method of fitting 
martingale measures to the marginal distributions. 

1 Introduction 

We consider the problem of finding real valued martingales fitting given 
marginal distributions and show that, by restricting to a certain class of 
strong Markov processes, it can be done in a unique way. It is furthermore 
shown that this is the unique continuous method of matching any specified 
marginals by martingales. 

The existence of martingales with specified marginals has been previously 
studied by many authors. In particular, Strassen [i'.i] showed in 1965 that 
if (/Un)neN is a sequence of probability measures on the real numbers which 
have constant mean and are increasing in the convex order, then there is a 
martingale such that the law of X„ is fin- The property that ^„ is 

increasing in the convex order simply means that Hn{f) is increasing in n for 
every increasing convex function /, and the necessity of this condition follows 
easily from Jensen's inequality. This result was extended by Kellerer [11] 
in 1972 to the case where the marginal distributions nt and the martingale 
Xt are indexed by time t in M_|_. It was also shown that X can always be 
chosen to be Markov. 

More recently, this problem has been investigated in the context of pric- 
ing financial derivatives, where knowledge of the prices of vanilla call and 
put options provides an implied distribution for the underlying asset price 
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at future times. For example, assuming zero interest rates (for simplicity) 
the local volatility model constructs the asset price process as a solution to 
the stochastic differential equation 



where S is a Brownian motion, S is the asset price and a{t, x) is the local 
volatility. Then, as is well known (see [ ] and [ ]), if C{t,x) is the price 
of a vanilla call with strike price x and maturity t, the implied probability 
density of St is d'^C/dx'^ and the local volatilities can be recovered from the 
following forward equation 



Alternative methods of matching the implied marginal distributions have 
been considered, such as jump-diffusions in [1], stochastic volatility in [2] 
and models based on Levy processes in [ ;]. Also, [16] gives several construc- 
tions, including Skorokhod embedding and time-changed Brownian motion 
methods. 

In this paper, we provide a general way of matching marginal distribu- 
tions under very mild constraints. Other than the necessary conditions of 
having constant mean and being increasing in the convex order, the only 
further constraint placed on the marginals is that they be weakly contin- 
uous. That is, if t„ — > f then Htnif) — > /"*(/) for every continuous and 
bounded function /. It is shown that such marginals can be fitted in a 
unique way by a certain class of strong Markov martingales. As this class 
includes all martingale diffusions, the solution will coincide with the local 
volatility model when it applies. However, for marginals which are either 
not smooth or don't have strictly positive densities, different types of solu- 
tions are obtained which cannot be described by an S.D.E. such as (1.1). 
For example, jump processes and singular diffusions, as described in Section 



We also show that that the resulting map from the sets of marginals to 
the martingales is continuous. So, a small change to marginal distributions 
results in only a small change to the martingale measure matching these 
marginals. 

Furthermore, it is shown in Theorem 1.5 that not only is our method 
of fitting the marginals continuous, but it is the only possible continuous 
method. Consequently, any alternative approach (e.g., those described by 
[1], [2] and [3]) must either fail to fit, or very closely approximate, certain 
marginal distributions, or small changes in the marginals would lead to big 
changes in the resulting martingale measure. 

Let us now define the types of processes to be considered, which should 
include all continuous and strong Markov processes. However, there are 



(1.1) 



dSt = St(T{t, St) dBt 
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some marginal distributions which cannot be matched by any continuous 
process. For example, if P(0 < < 1) = for all times t and ^ {Xt < 0) 
decreases in t, then there must be a positive probability that X jumps from 
below to above 1. For this reason, we relax the continuity condition to 
obtain the following class of processes. 

Definition 1.1. Let X be a real valued stochastic process. Then, 

1. X is strong Markov if for every bounded, measurable g : M ^ M and 
every t £ M+ there exists a measurable f : M+ x M — > M such that 

f{T,Xr)=E[g{Xr+t)\J'r] 

for every finite stopping time r. 

2. X is almost-continuous if it is cddldg, continuous in probability and 
given any two independent cddldg processes Y, Z each with the same 
distribution as X and for every s < t £ M+ we have 

F {Ys < Zs, Yt > Zt and y„ / Z„ for every u G (s, t)) = 0. 

3. X is an almost-continuous diffusion if it is strong Markov and almost- 
continuous. 

In [13] it was shown that almost-continuous diffusions arise when taking 
limits of continuous diffusions in the sense of finite-dimensional distributions. 
Note that condition 2 is equivalent to saying that Y — Z cannot change sign 
without passing through zero, which is clearly true for continuous processes 
by the intermediate value theorem. In what follows, we often abbreviate 
'almost-continuous diffusion' to ACD. 

An alternative way of representing the marginal distributions fit which 
we make use of is through the function C{t,x) = J{y — x)j^ d^t{y)- The 
property that fit is increasing in the convex order is then equivalent to C{t, x) 
being an increasing function of t. Furthermore, the distribution functions 
are easily recovered from oo,x]) = 1 -|- (72(t, This leads to the 

following space of functions. 

Definition 1.2. Let CP be the set of functions C : M+ x M ^ M such that 

1. C{t,x) is convex in x and continuous and increasing in t. 

2. C{t, x) ^ as X ^ oo, for every t £ M+. 

3. There exists a real number a such that C{t,x) + x ^ a as x ^ — oo 
for every t £ M_|_. 
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Continuity of C{t, x) in t is just requiring the marginals fit to be weakly 
continuous in t, and the third condition is equivalent to them having a 
constant mean. The property that a process X has marginals consistent 
with some C £ CP can be expressed as 

(1.2) C{t,x)=E[{Xt-x)+] 

and, conversely, if X is a martingale which is continuous in probability 
then C given by (1.2) will be in the space CP. The notation used here is 
borrowed from the financial interpretation where C(t, x) are call prices, with 
maturity t and strike x, although we just make use of C G CP as convenient 
representations of martingale marginals both in the statements of the main 
results below and in the proofs later. 

We use the space of cadlag real valued processes (Skorokhod space) with 
coordinate process X on which to represent martingale measures. 

D = {cadlag functions uj : M_|- M} , 
X:R+xD^R, {t,uj) ^ Xt{u}) =io{t), 
J' = a{Xt:teR+), 
J^t = a{Xs:se[0,t]). 

Then, (D,.F) is a measurable space and X is a cadlag process adapted to 
the filtration J^f The existence and uniqueness of the martingale measure 
fitting given marginals is now stated. 

Theorem 1.3. For any C £ CP there exists a unique measure P on (D,jr) 
under which X is an ACD martingale and (1.2) is satisfied. 

See sections 3 and 4 for the proof of this result, which involves a weak 
compactness argument to construct the measure and applies a result from 
[13] concerning limits of almost-continuous diffusions. Then, a backward 
equation developed in [12] is applied to show uniqueness. 

Given any C G CP, the notation ¥c will be used for the unique ACD 
martingale measure matching the marginal distributions given by C. This 
defines a map C i— > ¥c, which we shall show is continuous under the appro- 
priate topologies. 

Let M(D) be the set of probability measures on (D,J^). A sequence 
(P„)„gi^ in M(D) converges to P in the sense of finite-dimensional distri- 
butions if and only if Ep^ [Z] — > Ep [Z] for every random variable Z of the 
form 

(1.3) Z = f{Xt,,...,XtJ 

for ti, . . . ,tm G M+ and continuous bounded / : M. 

We use the topology of pointwise convergence on CP, so Cn ^ C if and 
only if Cn{t,x) C{t,x) for all {t,x) G M+ x M. Note that this is slightly 
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stronger than weak convergence of the marginal distributions, which would 
be equivalent to convergence of C„(t, x) — Cn{t, y) to C{t, x) — C{t, y). The 
continuity result for the map from the marginals to the martingale measures 
is as follows. 

Theorem 1.4. For every C G CP denote the unique ACD martingale mea- 
sure given by Theorem 1.3 by Fq- Then, the function 

is continuous, under pointwise convergence on CP and convergence in the 
sense of finite-dimensional distributions on A^(D). 

So, given any sequence C„ E CP converging pointwise to C € CP then 
Epj^^ [Z] — > [Z] for every random variable Z of the form (1.3). The 
proof of this is left until Section 4. 

Not only is the ACD martingale measure fitting the marginal distribu- 
tions uniquely defined, but it is also the only way of fitting the marginals 
in a continuous way, as the following result states. Here, we again use the 
topology of pointwise convergence on CP and convergence in the sense of 
finite-dimensional distributions on A^(D). 

Theorem 1.5. Suppose that we have a continuous map from a dense subset 
S of CP to the martingale measures 

such that for every C £ S the equality Eq^ [{Xt — x)+] = C{t,x) is satisfied. 
Then Qc = ^c- 

In particular this shows that the choice of the class of almost-continuous 
diffusions used to fit the marginals is not arbitrary, but was in fact forced 
upon us. The proof of Theorem 1.5 is left until Section 5, where the idea 
is that there are certain marginal distributions for which there is only one 
possible martingale measure. These correspond to extremal elements of CP, 
and form a dense subset. 

We finally note that the fact that C{t,x) is continuous and monotonic 
in both t and x for every C G CP implies that pointwise convergence is the 
same as locally uniform convergence, and the topology is given by the metric 

(1.4) d{Ci,C2) =sup||Ci(t,x) -C72(t,x)| A2-l^|-* : {t,x) G M+ x m} . 

So, we are justified in only considering limits of sequences (rather than 
generalized sequences) in the explanations and proofs of theorems 1.4 and 
1.5. 
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2 Examples 

In this section we mention some examples to demonstrate the kinds of pro- 
cesses which can result from different properties of the marginal distribu- 
tions. 

2.1 Continuous diffusions 

If C{t,x) is strictly convex in x for every t > 0, then the support of Xt 
under the measure given by Theorem 1.3 will be all of M and, consequently, 
X will be a continuous process (see [ ] Lemma 1.4). If, furthermore, C is 
twice continuously differentiable then it can be shown that X is a solution 
to the stochastic differential equation 



for a Brownian motion B and with a given by the forward equation 



Then, if a{t,x) is Holder continuous of order 1/2, the Yamada-Watanabe 
theorem ([17] V, Theorem 40.1) says that (2.1) uniquely determines the law 
of X. This is the familiar situation covered by the local volatility model, 
and widely employed in finance (see [(>]). 

2.2 Jump processes 

Now suppose that the supports of the marginal distributions are contained 
in the set of integers. Then, the process X must be an integer valued pure 
jump process. Suppose furthermore that P {Xt = n) > for every t > 
and integer n. Then, the almost-continuous property says that X cannot 
jump past any integer values, so can only jump between successive integers. 
Therefore, X must be a piecewise constant process with jump sizes ±1. For 
example, it could be a symmetric Poisson process (i.e., the difference of two 
standard Poisson processes). 

More generally, jump processes can arise whenever the supports of the 
marginal distributions are not connected intervals. Consider, for example, 
a smooth C{t,x) in CP which is strictly convex in x, so that the conditions 
considered in Section 2.1 are satisfied. Then define C by 



The corresponding marginal distributions then assign zero probability to the 
interval (0,1). Under the resulting martingale measure the process X 
will behave like a continuous diffusion satisfying the SDE (2.1) whenever 
X > 1 OT X < 0. However, the points {0, 1} will act like reflecting barriers, 
compensated by X sometimes jumping across the interval (0, 1). 



(2.1) 




(2.2) 
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2.3 Singular diffusions 

Now suppose that C{t,x) is strictly convex in x for every t > so that, as 
in section 2.1, we can conclude that X is continuous under the associated 
ACD martingale measure. If, however, C{t,x) is not twice differentiable in 
X then the marginal distributions will not be continuous with respect to the 
Lebesgue measure, and X will not satisfy a stochastic differential equation 
such as (2.1). 

For example, suppose that C G CP satisfies all of the properties consid- 
ered in section 2.1 and define 

C{t,x) = ^C{t,x) + ^max(-x,0). 

Note that C{t,x) is not differentiable at x = and the corresponding 
marginal distributions have an atom at 0. Under the ACD martingale mea- 
sure P^, the process X will behave like a continuous diffusion satisfying (2.1) 
away from 0. However, ¥^{Xt = 0) = 1/2, so X is sticky at 0, spending a 
positive time there. If, furthermore, C(t, x) is strictly increasing in t then 
X will not be constant over any time intervals. 

Similarly, it is not difficult to construct marginal distributions so that 
Xt is rational with probability 1 and with support equal to M, resulting in 
continuous processes spending almost all their time in the rational numbers. 
This is the case with the Feller-McKean diffusion ([is] III. 23), and similar 
situations can arise as a limit of random walks with randomly generated 
rates (see [n]). 

3 Existence 

We show how ACD martingales can be constructed with specified marginal 
distributions by taking limits of processes which match the marginals at 
finite sets of times. A weak compactness argument is used to prove existence 
of the limit. 

For any set S Q ]R_|_, let M"^ consist of the real valued functions on 
S. We consider M'^ as a topological space using the topology of pointwise 
convergence, and denote its Borel u-algebra by . The weak topology on 
the probability measures on (R^ is the topology generated by the maps 
P I— > Ep [/] for all real valued continuous and bounded functions / on M"^. 
We denote the coordinate process on R'^ by X^ , 

X^ : S xR^ ^R, {t,Lj) ^ Xf (tj) = uj{t), 

which has natural filtration {J^f)teM.+ given by, 

jrf = (T (Xf : s G S, s < t) . 
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Then, for any measure P on (D,^) we use to denote the measure on 
(M'^,^'^) obtained from the law of Xt under P with t restricted to S. 

In particular, if S is countable then M'^ is a Polish space, as it has a 
countable dense subset consisting of those uj such that uj{t) is rational for 
all t G S" and zero for all but finitely many t, and the topology is given by a 
complete metric 

d{uj,uj') = ^2""'min(|u;(s„) - uj'{sn)\,l), 

n 

where S = {si, S2, ■ ■ ■}■ 

Furthermore, a sequence of probability measures P„ on (D, ^) converges 
to P in the sense of finite-dimensional distributions if and only if P^ —>■ 
weakly for every finite subset S of M+. 

We now prove the result that we need in order to be able to find limits of 
sequences of martingale measures. The idea here is to use weak compactness 
in order to pass to convergent subsequences. 

A set P of probability measures on a Polish space is said to be tight if 
for every e > there exists a compact set C with P (C) < e for all P G P, 
and P is then weakly compact. In particular, for any tight sequence of 
probability measures P^, there is a probability measure P and subsequence 
Fjik converging weakly to P (see [9] Theorem 15.39.) This allows us to find 
martingale measures with specified marginals as limits of sequences. 

Lemma 3.1. Let C £ CP and (P„)„gN be a sequence of martingale measures 
on (D, J^) such that Ep^ [(X^ — x)^] C{t,x). 

Then, there exists a subsequence P„j. and a martingale measure P on 
(D,.F) such that P„j. F in the sense of finite- dimensional distributions. 
Furthermore, X is a martingale under F, continuous in probability and sat- 
isfies Ep [{Xt-x)+] = C{t,x). 

Proof. First, choose any t G M-|- and e > 0. For every K > 

Fn {\Xt\ >K)< Ep„ [{Xt -K + l)+ + l + {Xt + K-l)+- {Xt + K)+] 
C{t, K-l) + l + C{t, 1-K)- C{t, -K). 

As this can be made arbitrarily small by making K large, we see that for 
every e > there exists a > such that P„ > K) < e for every n. 
Letting S = {si, S2, . . .} be a countable dense subset of M+ and e > 0, this 
shows that there exists a sequence Kn > such that 

Fn{\XsJ >Kn)<2-''e. 

Letting A be the compact set of all oj G M"^ satisfying |ci;(sn)| < Kn, 

oo oo 

P^(M^ \ A) < ^Fn{\XsJ >Kn)<Y, 2-"e = e. 

n=l n=l 
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So the sequence is tight and, by passing to a subsequence if necessary, 
we may assume that it convergence weakly to a probabihty measure Q on 

For every t £ S and x, y G M, weak convergence gives 
Eq [(Xf - x)+ - (Xf - y)+] = Um Ep„ [{Xt - x)+ - {Xt - y)+] 

= C{t,x)-C{t,y). 
Letting y increase to infinity and using dominated convergence, 
(3.1) Eq [(Xf -x)+] =C(t,a;). 

If s < t are in 5, Z : M"^ — > M is ^^-measurable, continuous and such that 
ZXg is bounded, and < Z < 1 then. 



Eq [Z(Xf - x)+] = hm Eps [Z(Xf - x)+] < Hm Eps [Z(Xf 



n— >oo 

< Hm Eps ((Xf - x)+ - (Xf - y)+)] + y) 
= Eq [Z ((Xf - x)+ - (Xf - y)+)] + C(t, y) 

Letting y increase to infinity and using dominated convergence shows that 
(xf — x)+ is a Q-submartingale. So, 



Eq [ZXf] = ^Hrn^IEQ ((Xf - x)+ + x)] 

< Hm Eq [Z ((Xf - x)+ + x)] = Eq [ZX, 



Therefore, X"^ is a Q-submartingale. Furthermore, as C G CP, (3.1) shows 
that Eq [xf ] is independent of t and X"^ is a Q-martingale. This ahows us 
to extend Xf to ah t E M+ using Xf = Eq [Xf |J^f] for any u > t in S. 
We now show that X'^ is continuous in probabihty. As it is a martingale, 
it has almost-sure left and right limits Xf_ , Xf,_ for every t G M+ (for t = 
set Xf_ = xf ). Assuming that G S", taking the difference of the right and 
left limits of equation (3.1) in t and using the continuity of C gives. 



= Eq [(xf+-x)+-(Xf_ 

Eo 1 |^S_ >x>Xf+ or Xf+ >x>Xf_. } I 



l{Xf_>x}(^t+ -^t-) + l{xS_>a;>Xf+ or Xf+>x>Xf_} 1^^*+ 



So, Q(XfL > X > Xf^) = Q(Xif_ > X > XfL) = for every x, showing that 
X^_ = Xf^. As X'^ is a martingale and right-continuous in probability, it 
has a cadlag version and, therefore, there is a measure P on (D, J'^) satisfying 
P"^ = Q. Furthermore, X is a martingale which is continuous in probability 
under P. Taking limits of i S 5 also shows that Ep [{Xt — x)+] = C{t, x). 
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It only remains to show that P„ ^ P in the sense of finite dimensional 
distributions. We use proof by contradiction, so suppose that this is not the 
case. Then there would exist a random variable Z of the form (1.3) for a 
finite subset F = {ti, . . . , tm.} of M+ for which Ep^ \Z\ does not converge to 
Ep \Z\ . Passing to a subsequence if necessary, we may suppose that 



for some e > and every n. Setting S' = S U F the above argument shows 
that, by passing to a further subsequence, there exists a measure ¥' on 
(D, JT) such that F^' (P')"^'- In particular, (P')'^ = hm„^oo Pf = IP"^ and, 
by right-continuity in t, it follows that P = P' and P^ — > P'^ contradicting 



Combining Lemma 3.1 with the results of [ >] gives the following, which 
will be used to construct ACD martingale measures with specified marginals 
by taking limits of measures matching the marginals at finitely many times. 

Lemma 3.2. Let P„ be a sequence of ACD martingale measures on (D,J-') 
and C G CP be such that Ep^ [{Xf — — > C{t,x). Then, there exists a 
subsequence P„j. converging in the sense of finite- dimensional distributions 
to an ACD martingale measure P satisfying Ep [{Xt — = C{t,x). 

Proof. First, Lemma 3.1 says that there exists a subsequence P„j. converging 
in the sense of finite-dimensional distributions to a martingale measure P 
satisfying (1.2) under which X is continuous in probability. However, Corol- 
lary 1.3 of [' ] states that under such a limit, X is an almost-continuous 
diffusion. □ 

To complete the proof of the existence of the ACD martingale measure, 
it just needs to be shown that it is possible to fit the marginals arbitrarily 
closely and Lemma 3.2 will provide us with the required limit. There are, 
however, many different ways in which we can go about this. For example, we 
could construct a diffusion as the solution of a stochastic differential equation 
to match smooth C G CP. Alternatively, finite state Markov chains could 
be used to approximate the marginals by finite distributions. However, one 
way to exactly match the marginals at any finite set of times is to use a 
Skorokhod embedding to time change a Brownian motion, as we describe 
now. This uses the methods described in ['"]■ 

Fix any C G CP and times to < ^i- Also, let /if be the associated 
one dimensional measures, satisfying J {y — x)^ dfitiu) = C{t,x). Then, 
define the distribution function Fi{x) = ^2(^1) ^+) + 1 = fJ-tiii—oo, x]). For 
u G (0, 1) set (3{u) = inf{x G M : ^1(2;) > n}, and let gu : [/3(n), 00) ^ M be 



(3.2) 



Ep„ [Z] > Ep [Z] + e 



(3.2). 



□ 



gu{x) = C{t^,p{u)) + (x - p{u)) (n - 1). 
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Figure 1: Fitting marginals to C{t,x) 



Then, a{u) > (3{u) is chosen to satisfy gu{a{u)) = C{tQ, a{u)). This 
uniquely defines a{u) when C{ti, P{u)) > C{to, P{u)) (see Figure 1), other- 
wise we set a{u) = (3{u). Also, set Fq-|^(x) = inf{u G (0,1) : a{u) > x}. 
This is right-continuous and increasing from to 1 , so is another distribution 
function. 

If 5 is a Brownian motion with initial distribution St = supg^^ Bg 

is its maximum process, then a stopping time r can be defined by 

(3.3) r = inf {t e M+ : Fi{Bt) < Fl^{St)] . 

Then B'^ is a uniformly integrable martingale and B^ has law equal to ■ 
See [Jij] for details (Proposition 2.2 and Corollary 2.1). 

Lemma 3.3. Let C E CP and Iq < ti €z M+. Then, there exists an ACD 
martingale X such that E [{Xt — x)^] = C(t, x) for t = tQ,ti. 

Proof. Let -B be a Brownian motion with initial measure fito and r be the 
stopping time (3.3). Let ^ : (to; ^i) — *■ 1^ be any continuous function increas- 
ing from — cxD to oo. For example, 9{t) = {ti — 1)~^ — (t — ^o)""^- We also set 
9(t) = —oo for t < to and 6{t) = oo for t > ti. Defining the stopping times 

Tt = mf{s eR+:Bs> e{t)} , 

the ACD martingale can then be constructed as Xt = BrtAr- The distri- 
butions of Xtg and Xt-^ are /ifg and uti ([10] Proposition 2.2). As B'^ is 
uniformly integrable, X will be a martingale. The paths of X are very 
simple — if T is the first time at which tt > t then 

r max{Bo,e{t)), if t < T, 

* \ Br, iit>T. 
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So X can only have a single jump at time T, at which X^^ > Xx- 

It only remains to show that X is an almost-continuous diffusion. Con- 
tinuity in probability is easy. As Xt- > Xt and E [Xt-] = E [Xt] by the 
martingale property, it follows that Xt- = Xt almost-surely. 

The stopping times A r are hitting times of the strong Markov process 
{Bt,St), where St = sup^^t^s, so the time changed process (XtjSrtAr) will 
also be strong Markov. However, St-^at = Xt if Xt > 9{t) and X is constant 
as soon as Xt < 0{t), so X must be strong Markov. 

We finally show that X is almost-continuous, so choose a cadlag process 
Y independent and identically distributed as X. Ys > Xg, and Yt < Xt 
for s < t then we can let T be the first time at which Yt < Xt, at which 
point we must have Yt- = Xt- = 0{T). In fact, T must be the first time at 
which Yt- = Xt-, and so is previsible. So, the martingale property gives 
E [Xt-] = E [Xt] and E [Yt-] = E [Yt] and, as Xt- > Xt, Yt- > Yt, we 
have Yt = Xt- □ 

This is easily extended to match the marginals at a finite set of times. 

Corollary 3.4. Let C G CP and A C M+ be finite. Then, there exists an 
ACD martingale measure P on such that Ep [(X^ — = C{t,x) 

for all t e A. 

Proof, li A = {t} is a single time, then P can be taken to be the measure 
under which Xt is independent of t with the required distribution. For 
A = {to < ti < • • • < t„} we use induction on n. Suppose that there is 
an ACD martingale measure Pi matching the required marginals at times 
tQ,ti, . . . ,tn-i. By Lemma 3.3 there is an ACD martingale measure P2 
matching the required marginals at times tn-i,tn- Noting that Xt^_^ has 
the same distribution under both Pi and P2, we can join these two measures 
together at time tn-i to get P, which is the unique measure on (D,.F) such 
that 

Ep [AB] = Ep, [^Ep, [B\Xt„_,]] = Ep, [Ep, [A\Xt„„,] B] 

for all bounded random variables A, B where A is -measurable and 
B is a{Xt : t > t„_i)-measurable. That X is an ACD martingale under 
P follows from the fact that it satisfies these properties over each of the 
intervals [0,t„_i] and [t„_i,oo). □ 

Finally for this section. Corollary 3.4 is applied to construct the required 
ACD martingale measure. 

Lemma 3.5. For every C £ CP there is an ACD martingale measure P on 
(D,.F) satisfying Ep [{Xt — x)^] = C{t,x). 

Proof. By Corollary 3.4, there is a sequence of ACD martingale measures 
P„ satisfying Ep^ [(X^/^ — = C{k/n,x) for k = 0,1,..., n. Then 
Ep„ [{Xt — x)^] — > C{t,x), so the existence of the ACD martingale mea- 
sure P follows from Lemma 3.2. □ 
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4 Uniqueness 

Using a generalized form of the backward Kolmogorov equation, it was 
shown in [12] that continuous and strong Markov martingales are uniquely 
determined by their marginal distributions. In this section we apply the 
arguments developed there to ACD martingales. The idea is that if X is 
a diffusion satisfying the stochastic differential equation (2.1) and it is as- 
sumed that C{t, x) = E [{Xf — x)-^] is twice continuously differentiable, then 
the backward Kolmogorov equation 

can be combined with the forward equation (2.2) to obtain the following 
martingale condition for f(t,Xt), 

df d^C dC d'^f ^ 

— 1 = 0. 

dt dx^ dt dx"^ 

This applies to twice continuously differentiable functions /. Multiplying 
by a smooth 6{t, x) with compact support in M>o x M, integration by parts 
gives 

dldCde _ mdldC _ dCdldf\ ^^^^ 
dx dx dt dx dx dt dx dx dt J 

This expression is defined for differentiable functions, and differentiability in 
t can be relaxed by replacing terms such as J - {df /dt) dt by the Lebesgue- 
Stieltjes integral / • dtf. The class of functions we consider is as follows. 

Definition 4.1. Denote by D the set of functions f : M+ x M — > M such that 

• /(t, x) is Lipschitz continuous in x and cddldg in t, 

• for every Kq < Ki € M and T € M-|- then JJ^^^Jq \dtf{t,x)\ dx < oo, 

• the left and right derivatives of f(t, x) with respect to x exist every- 
where. 

Also, let T>}<i be the functions f £ T> with compact support in M>o x M. 

In particular, if X is a cadlag martingale then C(t, x) = E [{Xt — x)-(_] 
will be convex in x and cadlag and increasing in t, so C E D. Using /~(t, x) 
to denote the left limit of /(t, x) in t, it was shown in [ ] (Lemma 2.2) that 
for f,g £ T> the partial derivatives df{t, x)/dx and df^{t, x)/dx exist almost 
everywhere with respect to the measure Jj ■ \dtg{t,x)\ dx. This enables the 
following definition to be made. Here, where we suppress the arguments of 
functions inside the integral signs they are understood to be (t, x). 
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Definition 4.2. For every f^g^T) define the linear map fJ-^^g] '■ "Dr ~ 



Comparing with the right hand side of (4.1), we expect the martingale 
condition to be = Oi which is indeed the case for continuous processes. 

However, to include discontinuous processes jump terms need to be added, 
leading to the following definition. 

Definition 4.3. Let X be a cddldg martingale. Then, for every f ^ T) 
define the linear map flj : 2?k ^ 



.t>o 

where C is defined by C{t, x) =K [{Xf — and 

(4.2) Ji'iej)^ r\f{t,x)-f{t,Xt) + {Xt-x)f;{t,x))e-{t,x)dx. 

JXt- 

The sum of J^^{0, f) will be integrable (see [12]), so flj is well defined. 
We have the following necessary martingale condition. 

Lemma 4.4. Let X be a cddldg martingale and f £ D. If f(t, Xt) is a 
martingale then fij = 0. 

See [12], Theorem 3.8. To prove Theorem 1.3, the converse of this state- 
ment needs to be shown and, for almost-continuous processes, we need to 
demonstrate that the jump terms Jf{9,f) can be eliminated to express the 
martingale condition solely in terms of Ai[/,c]- 

Fortunately, the jump terms will indeed drop out of the expression for 
flj as long as / is chosen such that f~{t, x) is linear in x across each of the 
connected components of the complement of the support of Xt. To show 
this, we make use of the marginal support, which was defined in [15] to be 
the set constructed from the supports of the marginal distributions of X as, 

MSupp(X) = {{t,x) G M+ X M : X G Supp(Xt)} . 

We now prove the following. 

Lemma 4.5. Let X be an ACD martingale and f ^T) satisfy f^^{t,x) = 
for every t > and x outside the support of Xt. Then, flj = 

Proof. It was shown in Corollary 4.8 of [15] that for an almost-continuous 
process, the set 



S = {{t, x) G M>o X M : Xs- < x < Xs oi Xg < x < Xs-} 
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is almost surely disjoint from the marginal support of X. Then, the condition 
of the lemma gives /~ = on 5", so 

(4.3) f-{t,x) - r{t,Xt) + f-{t,x){Xt -x) = 

for all (t, x) G S, almost surely. Finally, as there are only countably many 
times for which ^ f and X is continuous in probability, /(t, x) = f^{t, x) 
whenever Xt / Xt- and (4.3) gives Jt^{9, /) = 0. □ 

The following result shows that it will be enough to only consider func- 
tions f £ V satisfying /~ = outside the marginal support of X. 

Lemma 4.6. Let X be an ACD martingale, t > and g : M ^ M be convex 
with bounded derivative k < g' < K. Then, there exists an f : IR+ x M ^ M 
such that 

(4.4) f{s,Xs)=M[g{Xtm 

for every s < t, where f{s, x) is convex in x, right- continuous and decreasing 
in s, k < f^ < K and /^^ = outside the marginal support of X. 

Proof. First, by theorems 1.5 and 1.6 of [ ], f{s,x) can be chosen to be 
convex in x with derivative k < f^ < K and satisfying equation (4.4) for 
all s < t. Theorem 1.7 of [ : ] also says that f{s,x) is continuous on the 
marginal support of X over the range s < t. We can extend /(s, x) across any 
bounded connected component of ]R\Supp(Xs) by linear interpolation. Also, 
we can linearly extrapolate w.r.t. x above the supremum of Supp(X5) with 
gradient K and below the infimum with gradient k. This gives f,22{s, x) = 
outside the support of Xg. For simplicity, over s > t, fix f{s,x) < g{x) to 
be a linear function of x independent of s. Then, for s < u < t, Jensen's 
inequality gives 

f{s,Xs) =E[f{u,Xu)\Ts]> f{t,Xs), 

from which it follows that /(s, x) > f{u, x) for every x in the support of 
Xs- By linear interpolation and extrapolation of /(s, x) outside the support 
of Xs, this inequality is satisfied for all x. So, f{s,x) is decreasing in s. 
We now define /(s, x) to be the limit of /(n, x) as u strictly decreases to 
s. Then, / is convex in x, right-continuous and decreasing in s, and has 
derivative satisfying k < f^ ^ ^- By continuity, /(s, x) = f{s, x) whenever 
s < t and x is in the support of Xs, so / also satisfies (4.4). 

It only remains to show that f^^{s,x) = for s < t and x outside the 
support of Xs- First, suppose that x is in a bounded connected component 
(a, 6) of M \ Supp(Xs). As f{s,x) was chosen to be linearly interpolated 
outside the support of Xs, 

r{s, x) > f{s, x) = {b- {{b - x)f{s, a) + {x- a)f{s, b)) . 
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The reverse inequahty follows from the convexity of /~ and = / at points 
(s,a) and {s,b) (by continuity). So, = 0. 

Now suppose that x is in a connected component (a, oo) of M\Supp(Xs). 
As X is a martingale with Xg < a, we must have Xu < a for all u < s. But 
/(n, x) was chosen to be linearly extrapolated with gradient K over x > a, 



f {s,x) 



lim f{u, x) 



f-is,a) + K{x 



so we again have f^^{s,x) = 0. Finally, by the same argument, f^^{s,x) 
for x in a connected component {—oo,b) of M \ Supp(Xs). 



= 

□ 



For a general f £ D, f{t, Xt) need not be a semimartingale, so it 
will be necessary to consider the more general class of Dirichlet processes. 
These processes were introduced by Follmer in [7], generalized to the non- 
continuous case in [20] and [4], then further extended to the form required 
here in [I ] and [ ]. So using the notation of [ ], a process V is said to 
have zero continuous quadratic variation if its quadratic variation exists and 
satisfies [V]t = J2s<t^^s- For brevity, we say that a process is a z.c.q.v. 
process if it is cadlag, adapted and has zero continuous quadratic variation. 
A Dirichlet process is then defined as the sum of a semimartingale and a 
z.c.q.v. process. Quadratic covariations of Dirichlet processes are well de- 
fined, and for any semimartingale X and Dirichlet process Y the integral 
/ X- dY can be defined using integration by parts. 



Xs- dYs = XtYt - XqYq 



f 

Jo 



Ys-dXs-[X, Y]t. 



The following result then says that flj gives a measure of the drift of the 
Dirichlet process /(t, Xt) for any f £ D. Here, 'D(X) denotes the functions 
6 £ T> which decompose as 6{t, Xt) = Mt + Vt for a cadlag martingale M 
and finite variation process V satisfying 



E 



/ '^{\Xs-.\<K}\dVs\ 

Jo 



< OO 



for all t,K >0. 



Lemma 4.7. Let X be a cadlag martingale and f £ V. Then there is a 
unique decomposition 



(4.5) 



fit,Xt)=Mt + At 



for martingale M and previsible z.c.q.v. process A with Aq = 0. Further- 
more, for any £ I?k "with 9 £ 'D(X), 



E 


sup 


I 




_t>0 





9-{s,Xs-)dA, 



< OO 
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and, 



E 



(s,X,)dAs 



For the proof of this, see [12] lemmas 5.4 and 8.1. The proof of the 
martingale condition jlj = for ACD martingales will make use of Lemma 
4.7 to show that the process A in decomposition (4.5) is zero. This will be 
done by showing that the time reversed conditional variation, VJ^(A) below, 
is zero. First, the reverse filtration is 

= a {{Xs :s>t}U{AeJ':F{A)=0}), 

and Vx{A) is defined by 



VJciA) =supE 



Y.^k{A 



tk 



.k=l 



where the supremum is taken over all sequences < to < ti < • • • < in in 
and all ^^-measurable random variables \Zk\ < 1. The following result 
from [ - ] (Lemma 7.2) will be required to bound the variation of A. 

Lemma 4.8. Let X be a cddldg real valued process and A he a z.c.q.v. 

process such that supj>o \At\ is integrable and Vx{A) < oo. 

Suppose furthermore that At — Ag is -measurable for all t > s and 
that there exists a measurable u : M+ x M ^ M such that AAt = u(t, Xt) for 
all t > 0. Then, A has integrable variation satisfying 



E 



\dA\ 



< VRA). 



In order to apply Lemma 4.8 it will be necessary to calculate AA, which 
we do by making use of the quasi-left-continuity of X. Recall that a cadlag 
process X is quasi-left-continuous if Xr- = Xt for every previsible stopping 
time r. 



Lemma 4.9. Every ACD martingale is quasi-left- continuous. 

Proof. Let X be an ACD martingale, a > 0, and (7 : M ^ M be bounded 
and Lipschitz continuous with coefficient K. Then, there exists a bounded 
fa '■ K>o X M ^ M such that fa{t,x) is Lipschitz continuous in x with 
coefficient K and, 

(4.6) fa{t,Xt)=E[g{Xt+a)\H 

for every t > (see [15] Theorem 1.5 or [13] Lemma 4.3). By linearity, this 
extends to all stopping times taking only finitely many values in (0, 00). 
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If n < t and /i : M ^ M is continuous and bounded, then the continuity 
of X at t together with bounded convergence gives 

E[h{Xu)fa{t,Xt)] =E[h{Xu)9{Xt+a)] = limE[h{Xu)9{Xs+a)] 

= limE [h{Xu)fa{s, Xs)] = limE[h{Xu) fa{s,Xt)] . 

s—*t s— >t 

If we let u increase to t then E — will go to 0, and we see that 

the above equality also holds for u = t. As fa{t,x) is uniformly continuous 
in X, this shows that is continuous in t for every x in the support 

of Xt- So, fa is continuous on the marginal support of X. By lemmas 4.3 
and 4.4 of [15], the continuity in probability of X implies that the paths of 
{t,Xt) and {t,Xt-) lie in the marginal support of X at all times. So, by 
taking limits of stopping times which take finitely many values in (0,oo), 
equation (4.6) extends to all stopping times. Similarly, taking increasing 
limits of stopping times gives 

UT,X^^)=E[g{Xr+a-)\:Fr-] 

for all previsible stopping times r > 0. We now note that fait, Xt) di^t) 
in probability as a ^ 0. So, by uniform continuity in x, fa{t,x) g{x) on 
the marginal support of X, giving 

g{Xr-) = lim fa{r,Xr-) = limE[g{Xr+a-)\J'r-] =E[g{Xr)\J'r-] ■ 

So, if h{x) is any bounded measurable function and u{x,y) = h(x)g(y), we 
can multiply by h^Xr-) and take expectations, 

E[u{Xr-,Xr^)]=E[u{Xr^,Xr)]. 

Finally, by the monotone class lemma, this extends to all bounded measur- 
able functions u : ^ R and taking u{x, y) = l{x=y} gives P {Xr- = X^) = 
1. □ 

Finally, we prove the martingale condition fJ-u^c] — 0- 

Theorem 4.10. Let X be an ACD martingale and f £ satisfy /^(i, x) = 
for t > and x outside the support of Xf. Then, f{t,Xt) is a martingale 
if and only if fJ^ij^c] = 0- 

Proof. First, Lemma 4.5 gives jlj- = lJ-[f^c]- Then, if f{t,Xt) is a martingale 
Lemma 4.4 gives Ai[/,c'] = 0- Conversely, suppose that /i[/,c] = ^^'^ l^t 
f{t,Xt) = Mt + At be decomposition (4.5). Also, choose any twice con- 
tinuously differentiable 9 : R>o x R ^ R with compact support and set 
B = J 6{s,Xs~) dAs, which is a previsible z.c.q.v. process (see [12] Lemma 
5.5). We will show that the necessary conditions to apply Lemma 4.8 to B 
are satisfied. 
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First, E [AMt|.7v_] = for every previsible stopping time r, and the 
quasi-left-continuity of X gives 

^A^ = E [/(r,XO - r{r,Xr-)\rr-] = fir^X^-) - r(r,X,_). 

By previsible section (e.g., [ ] Corollary 4.11), it follows that this equality 
holds simultaneously at all times. In particular 

(4.7) ^Bt = d{t,Xt) {f{t,Xt) - r{t,Xt)) 

for all times at which AXt = 0. As / = /~ at all but countably many 
times, and X is continuous in probability, the right hand side of (4.7) equals 
when AX ^ 0. Also, as B is previsible and X is quasi-left-continuous, the 
left hand side is also zero when AXt 7^ 0. So, (4.7) holds for all times, and 
ABt is a function of (t,Xt), as required by Lemma 4.8. 

We now show that VJ^{B) = 0. To do this, first choose any T > and 
let 5f : M — > M be bounded and a difference of convex Lipschitz continuous 
functions. Lemma 4.6 says that there is an h £ D such that h~^{t,x) = 
for X outside the support of Xt and, 

h{t,Xt)=E[g{XTm 
for t < T. Then, integration by parts gives 

(4.8) h{t,Xt){Bt-Bs) = [ h-{u,X^_)e{u,Xu-)dA+ [ Audh{u,X^_) 

J s J s 

(see [12] equation 5.5 and Lemma 5.6). As h{t, Xt) is a martingale over t < T 
it is clear that h G T^{X). Also, 9 is twice continuously differentiable and 
Ito's formula shows that it is in "D^X) (see [12] Lemma 4.3). So, by Lemma 
8.2 of [12], the product Of is also in D^X). So, we can take expectations of 

(4.8) to get 

(4.9) E[g{XT){Bt - Bs)] =E[h{t,Xt){Bt - Bs)] = fiif,C]{l[s,t)h9) = 

for s < t < T. Here, we have used that fact that the final term on the right 
hand side of (4.8) is a local martingale and applied Lemma 4.7. Taking 
limits as T decreases to t, (4.9) holds for t = T. Then, the monotone class 
lemma shows that (4.9) holds for all bounded and measurable functions g. 
So, if Z is a bounded ^^-measurable random variable, the Markov property 
for X gives 

E [Zt{Bt - Bs)] = E [E [Z\Xt] {Bt - B^)] = 0, 
and, V^{B) = 0. 

Also, At — As (and therefore Bt — Bs) is ^^-measurable for all s < t 
([12] Corollary 7.5). So Lemma 4.8 shows that B has zero variation and is 
therefore constant. Then, as 9 is arbitrary, A must also be constant and 
f{t, Xt) = Mt is a martingale. □ 
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We end this section with the proofs of theorems 1.3 and 1.4. 

Proof of Theorem 1.3. Lemma 3.5 shows existence of the ACD martingale 
measure. To prove uniqueness, suppose that we have two such measures P 
and Q. Choose any t > and convex and Lipschitz continuous 5^ : M ^ M. 
By Lemma 4.6 there is an / G P such that /~ = outside the marginal 
support of X and satisfying 

=Ep [g{Xt)\rs] 

for all s < t. Letting s increase to t gives f~(t,Xt) = g{Xt) P-a.s.. So, we 
may set /(s, x) = x) for all s > t, and /(s, X*) will be a P-martingale. 

However, it is clear that the stopped process X* is also an ACD mar- 
tingale under both P and Q. So, setting C{s, x) = C{s A t, x) we can apply 
Theorem 4.10 to get fJ'[f c] = 0. Another application of Theorem 4.10 shows 
that f{s, Xl) is also a Q-martingale and, 

f{s,Xs) = Eq [f{t,Xtm = Eq [g{Xtm ■ 

So X has the same pairwise distributions under both P and Q. As they are 
Markov measures with common initial distribution, this gives P = Q. □ 

Proof of Theorem 1.4- Choose any C G CP and sequence Cn G CP such 
that Cn — > C. We use proof by contradiction to show that Pc„ Pc; so 
suppose that this is false. Then there would exist an e > and a random 
variable Z of the form (1.3) such that 

(4.10) Ep^ JZ] > Ep^ [Z] + e 

infinitely often. By passing to a subsequence if necessary, we may suppose 
that this inequality holds for every n. Then, Lemma 3.2 says that by passing 
to a further subsequence we have Pc^ P for an ACD martingale measure 
P satisfying (1.2). So, the uniqueness part of Theorem 1.3 gives P = Pc, 
and Pc„ Pc contradicting (4.10). □ 

5 Extremal Marginals 

Note that the space CP is a convex subset of the space of all real-valued 
functions on M-|- x M. For any convex subset of a vector space there is a 
concept of extremal points — they are the points which cannot be expressed 
as a convex combination of other elements in the set. More precisely, for a 
convex subset 5 of a vector space V, an element x £ S is said to be extremal 
if given any y,z £ S and any A G (0, 1) such that x = Xy + {1 — X)z then 
y = z = x. We shall call an element of CP extremal if it is extremal among 
the convex set of all other elements of CP with the same values at time 0. 



MARTINGALE MARGINALS 



21 



Definition 5.1. A C ^ CP is extremal if given any Ci,C2 S CP and 
A G (0,1) such that Ci{0,x) = C2{0,x) = C{0,x) (Vx G M) and C = 
XCi + (1 - X)C2 then Ci = C2 = C. 

It can then be shown that there is a unique martingale measure fitting 
the marginals given by an extremal element of CP. 

Lemma 5.2. Let C G CP he extremal. Then, there exists a unique martin- 
gale measure P on (D, !F) satisfying Ep [[Xt — x)j^] = C(t, x). 

Proof. As Theorem 1.3 says that such a martingale measure exists, we only 
need to prove uniqueness. Let us start by showing that any such measure 
P is Markov. Given any T > and .Fj-measurable random variable Z with 
< Z < 1, define Ci, C2 by Ci(t, x) = C2{t, x) = C{t, x) for t < T and, 

(5.1) Ci{t, x) = Ep [ZEp [{Xt - x)+\Xt]] + Ep [(1 - Z){Xt - x)+] , 
C2{t, x) = Ep [(1 - Z)Ep [{Xt - x) + \Xt]] + Ep [Z{Xt - x)+] 

for t>T. It is easily checked that Ci,C2 are in CP, Ci(0,x) = 6*2(0, x) = 
C(0,x) and (Ci + C2)/2 = C As C is extremal, this imphes Ci = C. For 
t >T, putting Ci{t,x) = Ep [{Xt — x)^] into (5.1) gives 

Ep [Z{Xt -x)+]= Ep [ZEp [{Xt -x)+\Xt]], 

so Ep [{Xt -x)+\J't]= Ep [{Xt -x)+\Xt], and X is indeed Markov. 

Let us now suppose that P, Q are two such martingale measures. Choose 
any T > and X-p- measurable random variable Z with < Z < 1, and 
define Ci,C2 by Ci{t,x) = C2{t,x) = C{t,x) for t < T and, 

(5.2) Ci{t, x) = Ep [Z{Xt - x)+] + Eq [(1 - Z){Xt - x)+] , 
C2{t, x) = Eq [Z{Xt - x)+] + Ep [(1 - Z){Xt - x)+] 

for t >T. Again, it is easily checked that Ci, C2 are in CP with Ci(0, x) = 
6*2(0, x) = C{0,x) and (Ci + C2)/2 = C. As C is extremal this implies 
Ci = C. For t > T, putting Ci{t,x) = Eq [{Xt - x)+] into (5.2) gives 

Ep [Z{Xt - x)+] = Eq [Z{Xt - x)+] 

for all t >T. Therefore, P and Q are Markov measures for X with the same 
pairwise and initial distributions, so P = Q. □ 

Theorem 1.5 will follow once it is shown that the extremal elements are 
dense in CP. For C £ CP we will construct extremal elements of CP which 
match C any given increasing sequence of times. We start by showing that 
there is an extremal element matching C at two times. The method used 
here corresponds to the ACD martingales constructed in Section 3 (Lemma 
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3.3), which was based on the Skorokhod embedding described in [10], so let 
us recall some of the definitions from Section 3. 

Given a C G CP let nt be the corresponding marginal distributions 
satisfying J{y — x)j^d^t{y) = C(t, x). Fixing to < ^ii define the distribution 
function Fi[x) = C^^{ti^x+) + 1 = ^u^j ((— oo, x]) and for u G (0,1) set 
I3{u) = inf{x G M : Fi(x) > ti}. For x > /3(n) set 

gu{x) = C(ti, /3(n)) + (x - (5{u)) {u - 1) 

and define a{u) > (3{u) by gu{a{u)) = C{tQ, a{u)). This uniquely defines 
a{u) whenever C(ti,/3(n)) > C{tQ, f3{u)) (see Figure 1), otherwise we take 
(3{u) = a{u). 

We define C : M+ x M ^ M by setting C{t, x) equal to C{to,x) for t < to, 
C(ti,x) for t > ti and, 

r c(ti,x), iix<p{u), 

(5.3) C{t,x) = < C(to,x), iix>a{u), 

[ gu{x), if /3{u) < X < a{u) 

for to < t < ti, where u is set to u{t) = (t — to)/(ti — to)- We show that this 
does indeed give an extremal element of CP. 

Lemma 5.3. Suppose that C G CP and to < ti G M_(.. Then C defined 
by (5.3) is an extremal element of CP such that C{t,x) equals C(to,x) for 
t < to and C(ti, x) for t > ti. 

Proof. For to < t < ti, C(t,x) is a convex function of x lying between 
C(to, x) and C(ti, x) (see Figure 1). To show that C G CP it just needs to be 
shown that it is continuous and increasing in t. So, pick any s < t G (to, ti). 
Then, u{s) < u{t) and /3(n(s)) < P{u{t)). For any x < f3{u{t)) this gives 

C{t,x) = C(ti,x) > C{s,x). 

Also, if X > [3{u{t)) then 

C(t,x) = max(5(„(i)(x),C(to,x)) > max(5(„(^)(x), C(to, x)) = C'(s,x). 

So, C{t,x) is increasing in t. 

Now choose any x, y G M with C(to, x) < y < C(ti, x), and choose 6 < x 
to minimize u = (C(ti,6) — y)/{x — b). To see that this exists, note that 
choosing b small enough so that C(ti, 6) — C{to,b) < y — C(to, x) gives 

(C(ti, 6) - y)/(x - 6) < (C(to, b) - C(to, x))/(x - 6) < 1 

but the limit as 6 ^ — oo is 1. So, by continuity, it must have a minimum. 

Choosing t such that u{t) = 1 — (C(ti, x) —y)/{b — x) gives C^^{ti., 6+) > 
u-l> C^{ti,b-), and it follows that C{t,x) = C(ti, 6) + (x- 6)(1 - u) = y. 
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Therefore, 1 1— C{t,x) maps the interval [to,ti] onto [C{to,x),C{ti,x)] and 
must be continuous. 

It only remains to show that C is extremal, so suppose that C = ACi + 
(1 - A)C2 for Ci, C2 G CP, A G (0, 1), and 

Ci{0,x) = 02(0, x) = 0(0, x) = CXto,x). 

In order to show that Ci = C2 = C we shall make repeated use of the 
simple fact that if a non-trivial convex combination of two increasing func- 
tions is constant, then those functions must also be constant. In particular, 
C{t,x) = C{ti,x) for all t > ti so we must also have Ci{t,x) = Ci{ti,x) for 
i = 1,2 and t > ti. Similarly, Ci{t,x) = Ci{0,x) = C{t(),x) for t < Iq. 

Now choose any t G (to,ii) and set (3 = (3{u{t)), a = a{u{t)). Then for 
X < P, the definition of C gives C{t,x) = C{ti,x), so it is also true that 
a(t,x) = C,(ti,,T) fori = 1,2. 

Also, if X > a then C{t,x) = C{tQ,x). Therefore Ci{t,x) = Ci{to,x) = 
C{to,x). Furthermore, for all x in a) 

X{CiUt,x+) + (1 - X){C2Ut,x+) = C,{t,x+) = u{t) - 1. 

As {Ci\^{t, x+) are increasing functions of x, this shows that they are con- 
stant as X runs through this interval. Therefore, Ci{t,x) are linear functions 
of X over a). So, we have shown that 

r C,(ti,x), ifx</3, 
Ci{t,x) = < C{to,x), if X > a, 

( {{a - x)Ci{h,0) + {x- P)C{to, a)) /{a - if /? < x < a. 

Now suppose that there exists an x G M such that Ci(ti,x) < C{ti,x). 
Choose t such that u{t) = C^^{t\,x+) + 1 and set (5 = f3{u{t)), a = a{u{t)). 
It follows that /3 < X < CK, so the fact that Ci(t,x) and C(t, x) are convex 
in X and increasing in t gives 

(Ci),,(ti,x+) > {C{t^,a)-Ci{ti,(5))/{a-(5) 
> {C{to,a)-C{h,P))/{a-(3) 
>C,{ti,x-). 

Taking the right hand limits in x gives (Ci)^2(*ii ^+) — ^^2(^1!^+)- 

Similarly, if Ci(ti,x) > C(ti,x) then C2(ti,x) < C(ti,x) and the above 
argument gives (C2),2(ii, x+) > C2{ti,x+). So (Ci),2(ti, x+) < (72(*i>2;+). 

This shows that the function ^(x) = (C(ti,x) — Ci(ti,x))^ has a non- 
positive derivative everywhere and must be decreasing, so the limit g{x) — > 
as X —* ±00 implies that g is identically 0. So Ci(ti,x) = C{ti,x), from 
which it follows that Ci = C. □ 

The previous lemma allows us to construct an extremal element of CP 
matching C at an increasing sequence of times. 
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Corollary 5.4. Suppose that C G CP and < ti < ■ ■ ■ ] oo are in M+. 

Then there is an extremal C £ CP such that C{tk,x) = C{tk-,x) for each k. 

Proof. Without loss of generality, we assume that to = 0. By Lemma 5.3 
there exists extremal Ck € CP such that Ck{t,x) equals C(tfc-i,a;) for t < 
tk-i and equals C{tig,x) for t > t^ {k = 1,2,...). Define C by C{t,x) = 
Ck{t,x) for tk-i <t<tk. 

It is clear that C G CP. It just needs to be shown that it is extremal. 
So suppose that 

(5.4) C = ACi + (1 - A)C2 

for A e (0,1), Ci,C2 G CP, and Ci{0,x) = C2{0,x) = C(0,x). We use 
induction on k to show that Ci{t,x) = C2{t,x) = C{t,x) for t < t^. For 
k = this is a required condition, so suppose that k > 1 and that this 
holds for all t < tk-i- If 9it) = (t At^) Vt^-i then C{e{t),x) = Ck{t,x) is 
extremal, and (5.4) gives Ck{t,x) = Ci{t,x) = C2{t,x) for tk-i <t <tk. 
So, by induction, C = Ci = C2. □ 

Finally we complete the proof of Theorem 1.5, showing that C 1— > P^; 
given by Theorem 1.3 is the unique continuous map to the martingale mea- 
sures and matching all possible sets of marginals. 

Proof of Theorem 1 . 5. Choose any C £ S and Z be a random variable of 
the form (1.3). We just need to show that Eq^ [Z] = Ep^ [Z]. 

First, Corollary 5.4 gives a sequence (C„)neN of extremal elements of CP 
such that Cn{k/n,x) = C{k/n,x) for all /c G N and, consequently, Cn — > C. 

Then, as S is dense in CP, there is a sequence {Cn,m)mm such that 
d{Cn,m, Cn) ^ as m — > oo, where d is the metric on CP given by (1.4). 

Now, fixing any n G N, Lemma 3.1 says that, by passing to a subsequence 
if necessary, there exists a martingale measure P satisfying equation (1.2) 
and such that m ^ ^ ™ the sense of finite-dimensional distributions as 
m — > oo. However C„ is extremal so, by Lemma 5.2, P = Pc,i- Therefore, 
for every n we can choose an G N such that 



diCn,m„,Cn) < 2"", Eq^^^^^ [Z] - Ep^^ [Z] 



< 2- 



In particular, d{C,Cn,mn) ^ d{C,Cn) + 2~" so the continuity of C 
shows that Qc„ m„ tends to Qc in the sense of finite dimensional dii 
tions as n goes to infinity. Similarly, Fc„ tends to Pc giving. 



|Epo [Z] - Eq^ [Z]\ = hni Ep^^ [Z] - Eq^^^^^ [Z] 



< lim 2"" = 0. □ 

n— >oo 
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